Universality in tightly bound 3-boson systems 
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The effects of two distinct operations of the elements of the symmetry groups of a Hamiltonian 
on a quantum state might be equivalent in some specific zones of coordinate space. Making use 
of the matrix representations of the groups, the equivalence leads to a set of homogeneous linear 
equations imposing on the wave functions. When the matrix of the equations is non-degenerate, the 
wave functions will appear as nodal surfaces in these zones. Therefore, the equivalence leads to the 
existence of inherent nodal structure in the quantum states. In this paper, trapped 3-boson systems 
with different types of interactions are studied. The structures of the tightly bound eigenstates have 
been analyzed systematically. The emphasis is placed to demonstrate the universality arising from 
the common inherent nodal structures. 
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A. Introduction 

It is well known that universality exists in weakly 
bound 3-body system, where the Efimov states emerge, [l], 
[2l These states arise from the inherent mysteries of quan- 
tum mechanics, and they are governed by a universal 
law, not depend on the details of dynamics. In this pa- 
per, another kind of universality that exists in tightly 
bound quantum mechanic few-body systems is revealed 
and confirmed numerically. 

Usually, the Hamiltonian of identical particles is in- 
variant under the operations of a set of symmetry groups 
Ga, Gp, ■■ (including the permutation group). [j-Q Con- 
sequently, the eigenstates {\E'i(X)} are classified accord- 
ing to the representations of these groups, where X de- 
notes a set of coordinates and z is a serial number. Let 
ga be an element of Ga , 9/3 be that of Gp (each of them 
may be an element of the direct product of groups), and 
S denotes a special zone in the high-dimensional coor- 
dinate space. When X e S, the operations of ga and 
(7/3 might be equivalent so that ga'i'iiX) — g^5'i(X). 
For an example, when S is the zone of regular triangles 
(RT), ga is the rotation about the normal of the trian- 
gle by 27r/3, and gp is a cyclic permutation of particles, 
then ga and gp are equivalent in 5. Making use of the 
representations of groups, the equivalence leads to the 
establishment of a set of homogeneous linear equations 

J:AD?M - Df,M]'fAX) = in S, where D^ga) 
is the matrix element of the representation. When the 
matrix of this set of equations is non-degenerate, the set 
\I'i(X) must be zero in S. In this case, S becomes a 
prohibited zone and the wave function appears as an in- 
herent nodal surface (INS). [7] Since the matrixes of rep- 
resentations are absolutely irrelevant to dynamics, the 
appearance of the INS is universal disregarding the kind 
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of systems (nuclear, atomic, or molecular) and the details 
of dynamic parameters. 

The inherent nodal structures of 3-boson systems 
(without a trap) has been studied by the authors 
previously. [1, Isi] Accordingly, the quantum states can 
be naturally classified according to their inherent nodal 
structures. [sl-fllj Due to the progress of techniques, a 
few atoms can be tightly bound in a magnetic or op- 
tical trap. Thus the study of these systems is mean- 
ingful. The present paper extends the previous research 
of the authors and is dedicated to trapped 3-body sys- 
tems. A more efi'ective and simpler approach is proposed 
for the analysis. In particular, in addition to the RT, 
the symmetry constraint associated with the isoceles tri- 
angle (1ST) and the collinear geometry (COL) has also 
been taken into account. By studying the spectra, the 
root mean square radii, the particle spatial distribution, 
the weights of Q-components (Q is the projection of the 
total orbital angular momentum L on the third body- 
axis), and the shape-densities, systematic knowledge on 
the quantum states has been obtained. The emphasis is 
placed on demonstrating the universality by comparing 
different 3-boson systems. In addition to boson systems, 
a short discussion on fermion systems is also given at the 
end. 



B. Hamiltonian, spectra, and the eigen-states 



'Corresponding author: stsbcg@mail.sysu.edu.cn 



It is assumed that the three particles are confined 
by an isotropic harmonic trap ^muj^rf. Due to the 
trap numerous bound states exist, therefore a system- 
atic analysis can be made, hw and ^hjmu) are used as 
units of energy and length. Three types of interactions 
are assumed: Yk{t) = 10(2e-(''/i-428)' _ e-(r/2.i05)=)^ 
V^{r) = 1000e-3'--40/r6 (r > L2) and V^{r) = VA^.2) 
(r < 1.2), and Vc(r) = (r > 1) and Vb(r) = 15 (r < 1). 
The type Va has a short-ranged character and was pre- 
viously used in nuclear physics for the a-particles, Vb 
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FIG. 1: Spectra of 3-body systems with L = to 4 and par- 
ity n = ±1. Only the first-states are involved. For bosons 
three types of interactions Va, Vb, and Vc are considered. For 
fermions (the spatial wave functions are totally antisymmet- 
ric) only Va is considered. The levels have been shifted so 
that all the ground states have energy zero. The energy unit 
for the two cases with Va is fiu. In the case with Vb or Vc the 
energy unit has been redefined so that the energy of the 1"*" 
state has the same value as the one with Va. In every case the 
states are divided into three groups. For bosons, when differ- 
ent interactions are used the members and sequence in each 
group are the same (the only exception is that 4~ is higher 
than 4"^ when Vc is used). 



belongs to the Van der Waals type for atoms, while Vc 
is just a hard core potential. In fact, the interactions are 
chosen quite arbitrary, just to show the inherent univer- 
sality among different systems. The Hamiltonian is 
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(J = A,B, or C). 
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A set of basis functions is introduced to diagonalize the 
Hamiltonian to obtain the spectra and the eigenstates. 
The details are given in the appendix. Note that both Va 
and T^B contain a minimum, and therefore the low-lying 
states will pursue a better geometry so that the particles 
are appropriately separated and the total potential en- 
ergy can be thereby lower. This pursuit is not obvious 
for the case with Vc, where no minimum is contained. 

When the cm. motion is removed an eigenstate is 
denoted as where M is the Z-component of L, 11 

the parity, and i denotes the i-th state of a L'^-series. 
The 1 = 1 state (the lowest in a series) is called a first- 
state. The label i is usually ignored when i = 1. For the 
three types of interaction the spectra of nine first-states 
are plotted in Fig. [T] (the = 0~ state does not exist 
due to the "Rule 1" given below), where the spectrum 
at the right is for fermions and will be concerned later. 
It turns out that the three spectra at the left for bosons 
are similar. The levels can be similarly divided into three 
groups. In a group, the structures of the states are similar 
as shown below, and those with a larger L are higher due 



FIG. 2: The ratio r,n,s(i")/rons(0+), of the Lf states. 
Square is for Va, circle for Vb, and triangle for Ve- 



to the rotation energy. Many common features existing 
in the three spectra (say, 1+ is particularly high, 3+ is 
much higher than 3~, and 2" is much higher than 2+, 
etc.) can not be explained simply by dynamics. It implies 
the existence of a more fundamental reason. 

The root mean square radius r^msiL^) is ex- 
tracted from j^.f to evaluate the size. The ratios 
?'rms(-^'^)/?'rms(0^) are shown in Fig. [2] where the three 
dotted lines guiding the eyes go up and down in a syn- 
chronous way. For all the three types of interactions, the 
sizes of 0+, 2"*", and 3~ are relatively smaller, while the 
sizes of 1"^ is the largest. From dynamics, the 4+ and 
4~ states are expected to be larger in size because the 
particles are pushing out by a stronger centrifugal force, 
but in fact not. It demonstrates again the existence of a 
more fundamental reason. 

Let k' be a unit vector vertical to the plane of particles. 
For L states, the relative orientation of k' and L is 
a noticeable feature. Let a body frame S' with the base- 
vectors i' — j' — k' be introduced, where j' lies along 
R = ~ (ri + r"2)/2. Then, the eigenstates can be 
expressed as ^'"^.(X) = Eg D^j,,{n)^%{X') where X 
denotes the set of coordinates relative to a fixed frame S, 
and X' denotes the set relative to E', denotes the three 
Euler angles responsible for the transformation from E' 
to E.[1J] Starting from the normality 1 = (^l.a/|*l,m), 
we carry out first the integration over fl. Then we have 



1 = 



E 



2L+ 1 



dX'\^lQiX')\'^J2^Q, (2) 



where dX'^i?^drdi?sin6'd6', 6 is the angle between r = 
r2 — ri and R. Wq is the weight of the Q-component 
^J^g. Note that Q is the projection of L along k' . Thus, 
in ^2^, L is essentially lying along k' and therefore is 
essentially vertical to the plane of particles. Whereas in 
^2o' ^ is essentially lying on the i' — j' plane and there- 
fore is coplanar with the particles. Thus, the relative 
orientation between L and k' can be understood from 
Wq. 
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Carrying out the integration over dX' numerically, Wq 
can be obtained. Since Wq = W-q, we define Wq — 
2Wq (if Q ^ 0), or = PT/q (if Q = 0). They are given in 
Tab.Hl 



TABLE I: The weights of the |(5| -components obtained by 
using Va- 
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Wq are far from uniform but concentrated in one 
I component (except in 4+). When other interactions 
are used the qualitative features of Wq remain unchanged 
(say, for Va, Vb and Vc, respectively, 3~ has W3 = 0.984, 
0.991, and 0.972). Thus, again, there is a fundamental 
reason beyond dynamics. 



C. Symmetry constraint and the classification of 

states 



by the lower states. Among the nine states under con- 
sideration, the RT-ac components can be contained in 
the 0+, 2+, 4+, 3^ and 4^ states. They are called the 
RT-ac states and are dominated by the RT-ac compo- 
nent with IQI = (3) if n = 1 (-1) (as shown in Table 
1), thereby their energies can be much lower and they 
constitute the lowest group in the spectrum. Their wave 
functions are distributed surrounding a RT resulting in 
having a smaller size. The domination of the RT-ac com- 
ponent explains why Wq or W3 is particularly large in the 
RT-ac states. 

For l+j the |Q| = 1 component is prohibited by "Rule 
1" , while in the Q — component the 1ST is prohibited 
by " Rule 3" . The prohibition of all the 1ST causes serious 
consequence. It implies that the wave function is expelled 
from the important zones of coordinate space (see below) 
resulting in having a very high energy and a large size. 
This explains why this state belongs to the highest group 
in the spectra. 

For 1~ and 2^, only \Q\ = 1 is allowed due to "Rule 
1" . This component is IST-ac but RT-inaccessible (RT- 
inac). Accordingly, these states are classified to the sec- 
ond group. For 3"'", both |Q| = and 2 components are 
allowed. Due to "Rule 3", the former is IST-inac and 
therefore extremely unfavorable. Thus the weights are 
greatly concentrated in W2 as shown in Table 1. Since 
the IQI = 2 component is also IST-ac and RT-inac, 3+ 
belongs also to the second group. 

Since the classification is based on the accessibility of 
the RT and 1ST originating from symmetry constraint, 
it is universal for all kinds of 3-boson systems. 



To clarify the fundamental reason, we study the effect 
of equivalent operations. Let the operator of a rotation 
about an axis a by an angle /3 be denoted as 7?.^. When 
a is k' which is vertical to the plane of particles, the ro- 
tation by TT is equivalent to an inversion T with respect to 
the cm.. Thus we have 7^J'*"Q = e-^''^ = 
It leads to the "Rule 1: only Q even (odd) components are 
allowed for parity even (odd) states^^ as shown in the ta- 



ble. Furthermore, at a RT, T^2tt/3'^ lq 
^cyc^'iQ, where Pcyc denotes a cyclic permutation of par- 
ticles and will cause no effect. Thus, we have the "Rule 
2: ^2 Q is nonzero at a RT only if Q — 0, ±3, ±6, • • • ". 
When the three particles turn out to form an 1ST (in- 
cluding the RT as a special case) with the symmetric 
axis V lying on the i'-j' plane, we have the equivalence 
TZ'^'^^Q = Pij^^Q' where Py denotes the interchange 
of the two particles at the bottom of the 1ST. It leads 
to e-*25Q(„i)L+Q ^-n _Q = e*gg, where S is the angle 

between v and j' , and e = ±1 depends on the statistics. 
Thus, for bosons and for Q = 0, we have the "Rule 3: 
Q is nonzero at an 1ST only if L is even" . With the 
three rules, when L < 4, the RT- accessible components 
are ^-gzi^.^^o and^-nr-d^d ^g. 

Since all the three bonds can be optimized in a RT, 
the RT-accessible (RT-ac) components are much favored 



D. One-body densities 

Since the weights Wq are related to the relative ori- 
entation between fc' and L, these weights are related to 
the one-body densities pi which is an observable. Let us 
go back to the fixed-frame and deal with the eigenstates 
(r, R). From the normality, we have 



1 - / drdR\^%\ 



= J dr^piirs), 



(3) 



where = is the position vector of the third particle 
relative to the cm., and ipilr^) is the probability den- 
sity of finding a particle (disregarding which one) at r^. 
Let M — L, it implies that L is given along the Z-axis. 
The pi extracted from ^E"^^ is plotted in Fig. [3l For 3~, 
since this state is dominated by W3, k' is essentially ly- 
ing along L (i.e., along the Z-axis at the choice M = L). 
Therefore, it is expected that the pi of 3^ would be dis- 
tributed close to the X-Y plane. This presumption is 
confirmed in Fig. [31 For 4~, W3 is also dominant, thus 
the direction of L is close to but deviates a little from the 
Z-axis. Accordingly, the peak of pi of 4~ is close to but 
deviates from 6*3 = 90°. 3+ is dominated by W2 ■ Ac- 
cordingly, the peak is farther apart from 90°. Whereas 
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FIG. 3: One-body densities pi(7'3,6'3) of the ^f^i states cal- 
culated by using Va and with rs given at rrms- ^3 is the polar 
angle. P\{t^,6^) — pi(r3, vr — 6^3), and pi does not depend on 
the azimuthal angle. 

2+ is dominated by Wq. Therefore k' is essentially ver- 
tical to the Z-axis, thus the plane of the particles prefer 
to be coplanar with the Z-axis if \M\ — L, and the dis- 
tribution is no more close to the equator. In fact, the pi 
of 2+ is peaked at 63 — 0. 

Although Fig. [3] is from Va, the feature is common for 
various 3-boson systems. For examples, when Va, Vb, 
and Vc are used, the pi of 4~ is peaked at 63°, 63°, and 
66°, respectively, and the pi of 3+ is peaked at 50°, 49°, 
and 49°. 



E. Shape-densities 

Let us study the structures of the states in more 
detail in the body-frame. We define the hyper-radius 

f) = ■^/|rM-~|i?^. This quantity is invariant under the 
transformation between different sets of Jacobi coordi- 
nates. Therefore, it is suitable to d escrib e th e size. Fur- 
thermore, we define tan/3 = A/2/3i?/(\/l/2r). When r 
and R are replaced by f) and (3 as arguments, Eq. ^ can 
be rewritten as 

1 = J dt)dS PsUX'), (4) 

where d5 = (i+tan^ ^^^^^ ^^ d/S sm0d9 is an in- 
finitesimal deformation, and 

is the probability density that the system has a given size 
and a given shape (while an integration over the orienta- 
tion of the shape has already been carried out). Therefore 
Psha.{X') is called the shape-density which relates directly 
to the structure. 
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FIG. 4: Geometries in the ((/>, R/r) subspace. Where both the 

curves BAC and B' AC correspond to the 1ST. These curves 
fulfill the equation R/r — ^\/5 — 4F, and F = ^(cos^0± 

\/(3 -f- sin^ (j)) sin^ (f)). The vertical line (f> — Q corresponds 
also to the 1ST. The intersection A ((^ = 0, R/r = \/3/2) cor- 
responds just to a RT. Both the vertical lines (j) — ±90° corre- 
spond to coUinear structures COL. In particular, the point B 
and B' together with the horizontal line R/r = correspond 
to a symmetric COL, namely, a particle sits at the middle of 
the other two. The point C and C' correspond also to a COL 
but with two particles overlap with each other. 



Let (j) be the azimuthal angle of r lying on the i'-j' 
plane. Since R is given lying along j\ cj) — ■n/2 — 9. 
Let [) be given at vermis (Let ri be the distance of the 
«-th particle from the cm.. Then f) = \/rf + + , 
and the average 1} — (|rf + + rl\)^^'^ = VSri-ms)- Ex- 
amples of psha as a function of (j) and R/r = -^tan/3 
are given below. The geometry associated with the ar- 
guments {<j>,R/r) is shown in Fig.H) 

When Va is used, the psha of the RT-ac first-states 
are very similar with each other (except 4+) as shown in 
Fig. [5^ and [Sh. The psha of 4+ is shown in [St. All of 
them are peaked at {(j), R/r) — (0, \/3/2) associated with 
a RT as expected. In 5c the peak is extended along three 
directions each leading to a COL. All the three RT-inac 
and IST-ac first-states are similar as shown in [SJi, where 
a node appears at the RT and the three peaks are all 
associated with a flattened 1ST with particles 3, 1, and 
2 at the top, respectively. The top angles of the 1ST 
are 94°, 101°, and 105°, respectively for 3+, 1", and 2". 
The psha of the IST-inac state 1"*" is shown in[5f, where 
a number of nodal lines exist in the (0, R/r) subspace, 
and the wave functions are expelled and compressed into 
six separate regions. The appearance of so many nodal 
lines implies the existence of a very strong inherent os- 
cillation. This explains why 1"*" is very high in energy. 
From the different patterns of the psha, we know that the 
classification based on the accessibility of geometries is 
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valid. 

For a COL, let k' be set along the line of particles, then 
the rotation about k' by any angle causes nothing. Hence 
the COL exists only in the Q = component (it implies 
that L is essentially vertical to k'). Besides, a rotation 
about j' by TT is equivalent to an inversion. This leads 
to (-l)^vlin^ ^ n^P^p. Accordingly, only the states with 
(—1)^ = n are allowed to have the COL. Furthermore, 
for a symmetric COL, an inversion is equivalent to an 
interchange of the two particles at the ends. Thus, the 
symmetric COL exists only in positive parity states (for 
bosons). When a COL is accessible but the symmetric 
COL is not, the COL is not stable because an inherent 
node is contained and a linear oscillation is inherently 
excited. On the other hand, from pure dynamics, a COL 
will have a larger moment of inertia. Comparing with 
the RT, a COL will have a higher interaction energy but 
a smaller rotation energy. Hence, this structure will be 
preferred by the states with a larger L so that the rotation 
energy can be considerably reduced. With this analysis, 
among the nine states, obviously, the COL will be more 
preferred by 4+. This explains the extension found in 
[St, but not in [S^ and \5jp. It is expected that the COL 
will definitely appear in i > 2 states as shown in [S]; for 
4^, where the most probable shape is a symmetric COL. 
Incidentally, 1^ as shown in 5f is not only IST-inac but 
also COL-inac, this explains why so many nodal lines 
emerge. 

When Vb is used, the qualitative features of Fig. Ore- 
main unchanged. Note that all the first states will do 
their best to pursue a structure favorable in binding and 
free from the symmetry constraint. This leads to the 
similarity. When Vc is used, the total potential energy is 
not sensitive to the details of geometry because the inter- 
action does not contain a minimum. Hence, the pursuit 
to a better geometry is less anxious. This is shown in [6^ 
(Comparing with[5ji, the 1ST is no more well kept) and 
in|B}D (where the distribution extends towards the COL. 
This tendency is not explicit in 2'^ if Va or Vb is used). 



F. Fermion systems 

For fermion systems, we consider the case that the spa- 
tial wave function is totally antisymmetric against par- 
ticle permutation. Then the above Rule 1 and 2 hold 
exactly while Rule 3 becomes "VP^ q is nonzero at an 
1ST only if L is odd". With the three rules of fermions, 
when L < 4, the RT- ac states are 1+, 3+, 3~, and 4". 
The RT-inac but IST-ac states are 1", 2+, 2", and 4+. 
The IST-inac state is O"''. The spectrum from numerical 
calculation by using Va is shown at the right of Fig. [1] 
where the classification into three groups is revealed. As 
bosons, only the states with (—1)^ = H are allowed to 
access the COL. However, the symmetric COL is acces- 
sible only to H = —1 states. From an analysis of the 
Q-components, we know that 3~ has its RT lying on the 
i'-j' plane, while the RT of 3^ is upstanding. Further- 




FIG. 5: The shape-densities psha of some selected Li first- 
states (an exception is [S^ for the second-state) obtained by 
using Va- For each state the hyper- radius t) is given at VSr^ms, 
and the orientation of the shape has been integrated. The 
maxima of psha are marked by a cross and their values are 
given inside the parenthesis in each panel. The values of the 
contours form an arithmetic series and decrease to zero. 



more 3~ is symmetric-COL-ac, while 3"'" is not. This 
explains why 3~ is lower than 3^ even though they be- 
long to the same group. Similarly, the 1ST in 2+ is lying, 
that in 2 ~ is not (due to the prohibition of Q even com- 
ponent). Besides, 2+ is COL-ac while 2~ is not. This 
explains why 2"'' is lower than 2~ even though they be- 
long to the same group. 

Examples of psha of the 3-fermion states are shown 
in Fig. ^ and |6ji. Fig. |6}; is similar to [5f where the 
wave function is absolutely expelled from the lines of 
1ST. However, O"*" of fermions is COL-ac (but symmetric 
COL-inac) while 1"^ of bosons is not. Therefore, un- 
der the same interaction the 0+ of fermions is lower (see 
Fig. [!}, and its most probable shape is a non-symmetric 
COL with an inherent strong linear oscillation (due to 
the nodal lines) . Comparing Fig. ^ and \5jp (both for 
4''" but for fermions and bosons, respectively), the wave 
function of the former is expelled from the RT and also 
from the symmetric COL, while the latter is attracted by 
the COL. Thus the former is higher. 
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FIG. 6: psha obtained by using Vc for bosons (a and b), and 
Psha obtained by using Va for fermions (c and d). Refer to 
Fig. El 



G. Final remarks 



In conclusion, the features of the tightly bound 3-boson 
states have been studied systematically. The universal 
existence of inherent nodal structures has been demon- 
strated analytically and numerically. The existence is 
irrelevant to the strength of the trap. The quantum 
states belonging to the same set of representations of 
the symmetry groups are constrained in the same way 
disregarding the kind of particles. Since the constraint 
arising from symmetry is absolutely not violable, the ef- 
fect of symmetry overtakes dynamics and is decisive. In 
particular, there are prohibited zones (points, hues, ...) 
in the coordinate space not accessible to specific states. 
This leads to the establishment of a universal classifica- 
tion scheme based on the accessibility of RT, 1ST, and 
COL, and leads to the extensive similarity among differ- 
ent systems. Although this paper is essentially dedicated 
to boson systems, universality exists also in fermions sys- 
tems due to the decisive role of symmetry constraint. 
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Appendix 

I. DIAGONALIZATION OF THE 
HAMILTONIAN AND AN EVALUATION OF 
THE ACCURACY 

When the cm. coordinates have been removed, the 
Hamiltonian for internal motion Hint can be expressed 
by a set of Jacobi coordinates r — r2 — ri, and R = 
r-3 - (ri + r2)/2. Let 

/i(//,r) = -i-V2 + i^r2, (6) 

which is the Hamiltonian of harmonic oscillation. Then 
Hint = h{l/2, r) + h{2/3, R) + - (7) 

Let us introduce a variational parameter 7. The eigen- 
states of h{'j, r) are denoted as (Pni{r) with eigenenergy 
2n + Z + 3/2, where n and I are the radial and angular 
quantum numbers, respectively. From 4'nii''^) * 
basis functions for the 3-body system is defined as 

1>lnLA/(123) ^ ifZlir)<f>'^:£iR)),M, (8) 

where la and lb are coupled to the total orbital angular 
momentum L and its Z-component M, H = (— 
the parity, k denotes ria, la, n-h, and lb- These functions 
should be further (anti)symmetrized, thus we define 

KuLM = J2 *fe,nLAfblP2P3), (9) 
P 

where the right side is a summation over the permuta- 
tions. Making use of the Talmi-Moshinsky coefficients, 
p^MlSl each term at the right can be expanded in terms 
of*fe',nLM(123). Say, 

n.ULMim - $:^?,^1>Z, nLM(123), (10) 

k' 

where the Talmi-Moshinsky coefficients A^^, can be ob- 
tained by using the method given in the ref. (ISj. Note 
that the set {l>^ ulaA been orthonormahzed. 

Thus a standard procedure is further needed to transform 

hlm} *° ^ orthonormahzed set {^Jyilm}' where q is 

a serial number. Each $^ ulm '-^^^ expanded in terms 

of {*fe,nLA/(123)}. Finally, the set {^^^^ulm} is used to 
diagonalize Hmt- Due to the symmetrization the matrix 
element of any pair of interaction is equal to that of the 
pair 1 and 2. Thus only r and R are involved in the cal- 
culation. This leads to a great convenience. When 11 and 
L are given, a series of states will be obtained after the 
diagonalization. The lowest one of the series is called the 
first state. Its energy should be minimized by adjusting 
7- 
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Due to having a trap, the bound states are highly local- 
ized, the above basis functions are suitable and will lead 
to a good convergency. Let Nab — '2{na+nb)+la+lb- The 
number of basis functions is constrained by Nab < ^max, 
the latter is given. As an example for the bosons, the 



energy of the lowest 0+ state with the interaction Va 
would be 2.06561, 2.06558, and 2.06557 hui, respectively, 
if N„iax = 16, 18, and 20. Since the emphasis is placed 
at the qualitative aspect, the convergency appears to be 
satisfactory. 
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